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A model  of  the  axially  symmetric  vortex  has  been 
proposed  by  Nicholson  and  Johnson  (1978),  which  describes 
the  vortex  in  terms  of  its  vorticity  distribution  as  a 
function  of  radius.  The  vortex  model  is  a composite  of 
three  lateral  regimes  in  which  the  vorticity  distribution 
is  a piecewise  smooth  and  continuous  function  of  radius. 

This  distribution  along  with  a hypothesized  divergence 
distribution  is  shown  in  Figure  la,  with  the  concurrent 
tangential  and  radial  velocity  distributions  in  lb. 

Nicholson  and  Johnson  (1978)  present  data  from 
hurricanes,  tornadoes  and  waterspouts  to  demonstrate  that 
the  model  is  a good  and  reasonable  approximation  to  the 
tangential  wind  above  the  boundary  layer.  This  distribution 
satisfies  Laplace's  equation,  hence  the  name,  the  zero 
Laplacian  vortex.  Examples  of  the  fit  of  the  velocity 
to  the  data  are  given  below  in  Figures  2a  and  2b. 

Satisfaction  of  Laplace's  equation  by  the  vor- 
ticity distribution  suggests  that  the  vorticity  may  be 
part  of  an  harmonic  function  (the  imaginary  part)  and  that 
the  real  part  of  the  function,  the  divergence,  also  satisfies 
Laplace's  equation,  while  observing  matching  boundary 
conditions . 
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This  assumption  implies  that  a distribution  of 
divergence  accompanies  the  vorticity  distribution.  Assuming, 
as  a first  approximation,  that  the  atmosphere  is  horizontally 
incompressible,  then  from  continuity  the  vertical  velocity 
distribution  may  also  be  postulated. 


Laplace's  equation  is  a classic  elliptic  equation 
which  has  complex  characteristics.  The  axially  symmetric 
solution  (linearity  in  Unr)  is  not  the  only  solution  to 
Laplace's  equation.  The  asymmetric  solution  is  the  product 
of  an  exponential  and  periodic  function  such  that  the 
Laplacian  in  Cartesian  coordinates  is  given  by 
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is  a Fourier  amplitude  dependent  on  n, 
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The  trick  is  to  ascertain  the  relevance  of  the  solution  to 
the  phenomena  of  the  asymmetric  vortex.  This  will  be 
addressed  in  the  following  section. 


Spiral  Cloud  Bands 


It  has  long  been  established  that  the  spiral 
cloud  bands  in  hurricanes  occur  in  the  form  of  equiangular 
or  logarithmic  spirals.  That  is  to  say 


tan  a = 


= = c°nstant 


where  a is  the  angle  a log  spiral  makes  with  a ray  emanating 
from  the  center  of  a polar  coordinate  system. 

Since  a log  spiral  satisfies  (4)  then  it  should 
plot  in  the  semilog  coordinates,  9,  dnr  as  a straight  line. 
Interestingly  enough  when  multiple  cloud  bands  are  plotted 
in  this  coordinate  system  they  not  only  form  straight  lines, 
but  are  parallel  to  one  another  and  exhibit  a certain  amount 
of  periodicity.  Since  the  bands  increase  in  intensity  in 
a direction  orthogonal  to  the  periodicity  one  wonders 
whether  they  may  exhibit  behavior  which  would  satisfy  (2) 
and  (3). 
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In  Cartesian  coordinates  it  is  axiomatic  that 
rotation  of  coordinates  implies  a direct  transformation  of 
the  form  of  the  Laplacian  from  the  old  coordinate  axes  to 
the  new  ones  (Figure  3).  The  Laplacian  before  rotation  is 
expressed  as 


V2  s — _ + — _ 

2 3x2  9y2 


and  after  rotation  as 


v2  = -it 

2 Tx^ 


In  polar  coordinates  the  Laplacian  is  written 

v2  = I -ir  _ i + 4 

2 r 9r  9r  r7  Te7 

By  multiplying  the  radial  term  by  r/r,  recombining  and 
factoring  out  1/r2  the  Laplacian  may  be  rewritten  as 


2 r- 


9 2.nr‘ 


Laplace's  equation  may  dispense  with  the  1/r2  to  give 


r2 V2  a 


9&nr- 
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Thus,  any  linear  solution  in  8 and  £nr  satisfies  (9)  as 
well  as  the  product  of  an  orthogonal  set  of  exponential 
and  periodic  functions.  As  can  be  seen  in  (6)  a rotation 
of  the  coordinate  system  9 , £nr  by  angle  a in  Figure  4 
would  give  a new  expression  for  Laplace's  equation 


where  Sr  and  S0  are  the  new  coordinate  axes  illustrated 
in  Figure  4.  This  rotation  brings  the  cloud  bands  parallel 
to  SQ  and  periodic  in  Sr<  Since  both  SQ  and  Sr  satisfy 
(4)  they  are  log  spirals, and  the  new  coordinate  system 

s s 

e 9 and  e r are  log  spiral  coordinates.  S and  SQ  are  in 

r w 

units  of  radians.  They  vary  in  the  model  from  0 to  2ir, 
increasing  to  it  and  decreasing  again  to  zero. 

The  scalar  transformations  between  0,  Jlnr  and  S0, 

S are 
r 

S = Unr  cosa  + Qsina  (11a) 

r 

Sq  * 8cosot  - £nr  sin  a (lib) 

Thus  (2)  may  be  rewritten  as 

a * y(n)e  v cos(nSr),  etc.,  being  exponential 
in  SQ  and  periodic  and  linear  in  Sr. 
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3.  Spiral  Space 
3.1  Polar  Coordinates 

Spiral  space  is  represented  pictorially  in  Figure 

5 below.  The  spirals  are  defined  by  S and  SQ  being 

r u 

constants . 


Sr  = K = inr  cosa  + sina 


Sq  = K ' = cosa  - £nr  sina 

the  orthogonal  log  spirals  which  define  spiral  space  are 
given  by 


r = exp 


K - sina 


for  constant  S , i.e.,  the  coordinate  along  which  SQ  only 
r y 

changes  and 


cosa  - K' 


for  constant  Sq  , i.e.,  the  coordinate  along  which  only  Sr 

changes.  Thus  any  point  in  spiral  space  may  be  defined 

in  terms  of  values  of  the  log  spirals  e^r  and  e^® , or 

simply  S and  Sa . 
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3.2  Spherical  Coordinates 

Spiral  space  may  also  be  defined  in  geographic  or 
spherical  coordinates,  so  that  all  algorithms  applicable 
to  polar  coordinates  may  be  applied  to  spherical  ones  as 
well.  The  two  dimer  ional  Laplacian  in  spherical  coordinates 
$ , X is  given  by 


v2  „ 1 _3  sin<{>  _3  + 1 

2 r2sin<f>  d<p  3<f>  r7  sin2(J>  3A: 


rearranging  terms  gives 
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Rotating  the  coordinate  system  2,n  tan  ^ , A through  the 
angle  a produces 


V2  = — 2 

2 r sin2<}> 


where  S^  and  S,  are  log  spirals  on  a sphere  conforming  to 
the  equation 
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The  elliptic  solution 


>(S0  + i Sr) 


may  be  equally  applied  to 


.(  sx  + 1 se) 


where 


S,  = ln(tan  77)  cos  a + A sin  a 
(J)  ^ 


and 


= Xcos  a - £n(tan  ;|)sin  a 


The  spherical  space  is  defined  in  terms  of  spherical  log 
spirals  so  that 
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p = 2 tan  (exp 
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for  constant  and 
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S = v cos  a - In  tan  (-£)  sin  a 
v 2 


4.  Transformations  in  Complex  Space 


The  transformation  from  polar  to  spiral  coordinates 
involves  the  successive  transformations  of  shortening, 
rotation,  scaling  and  stretching. 


By  considering  complex  space 


Z = re 


ie 


and  its  conjugate 


Z = re 


-ie 


then  by  the  transformation  of  shortening  (operating  on 
the  conjugate  in  the  northern  hemisphere)  we  obtain 


l nz  = 2.nr  - i0 


Scaling  by  i and  n where  i = / -1  and  n is  an  interger 
produces 
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(12a) 


(12b) 


(13) 


ni  Unz  = n(0  + i inr)  (semilog  coordinates)  (14) 
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Rotation  by  the  versor,  e (equivalent  to  rotating  the 
semilog  coordinates),  produces 


nieict  SLnz  = n(9  ± £nr)ela  = n(S  + iS  ) 

y r 


Restretching  to  obtain  real  space  curvilinear  coordinates 
produces 


= exp  n(0  + i £nr)  elc 


= exp 


i(SQ  + iSr) 


The  Laplacian  in  complex  coordinates  is  given  by 


(16a) 


(16b) 


V2  = 4 -1  -4 
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so  that 
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and  16a  and  b satisfy  Laplace's  equation 


Assuming  the  divergence  and  vorticity  are  the 
real  and  imaginary  parts  of  un  harmonic  function,  i.e.,  that 
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they  satisfy  the  Cauchy  Riemann  equations  which  in  spiral 
coordinates  are 


95  _ 9 C,  9C  96 

9Sr  - 9S0  ’ 9Sr  - - 9Sq 


(19) 


then  we  may  write  the  asymmetric  solution  for  divergence 
and  vorticity  as 


6 = en^®cos  n S , c,  = - ens®sin  n S 

n r sn  r 


nS< 


where  this  solution  is  F where  F = 6 - i£  so  that 


(20) 


Re(F)  = en(S0)  cos  S 


(21) 


and 


Im( 


F)  = -en^S9^  sin 


In  this  case  the  divergence  field  leads  the 
vorticity  field  by  tt  / 2 . Thus  the  full  solution  for 
symmetric  and  asymmetric  distributions  is 

00 

6 = 6 +i  Jlnr  + / y (n)en^®^cos  n S 
o 1 o o r 

n=l 


(21a) 
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c = C0  + iQ£nr  - ^ Y0(n)en^Se)sin(n  Sr)  (21b) 

n=l 

where  r is  dimensionless,  being  understood  as  r/rQ  where 
rQ  = 1L.  5q,  £ , i and  yq  are  constants  of  integration 

derived  from  boundary  conditions. 

5.  Covariances  in  the  Divergence  and  Vorticity  Equations 

In  spiral  coordinates  the  divergence  and 
vorticity  may  be  written: 

6 = exp  (-2(Sr  cos  a - SQ  sin  a))  x 

g 

jg—  exp  (Sr  cosa  - Sg  sin  a)Ug  (22a) 

. r r 

+ -jg-  exp  (Sr  cosa  - SQ  sin  a)  Ug  1 


and 
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£ = exp(-2(S  cosa  - SQ  sin  a))  x 
r a 
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l9Sr 


exp  (S  cosa  - S sina)  U,, 
r u o , 


(22b) 


exp  (S  cosa  - S0  sin 


in  a)  Ug 
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Since 


6 = y (n)en^S®^cos  nS 

no  r 


Cn  = Yo(n)etl(S0'>sin  n Sr 


then  by  substitution  and  partial  integration  both  Uc  and 

r 

U„  have  as  part  of  their  form 
b0 


Un  / e 


ax  cos  pxdx 


°n  e 


ax  sin  pxdx 


which  integrate  to 


ax. 

e ( a cos  px 
a2  + p2 


clX 

e (a  sin  px 


sin  px. 


cos  px. 


a2  + p2 


(24) 


respectively.  Hence  both  Uc  and  U„  are  out  of  phase  by 

r se 

a non  multiple  of  tt/2  with  the  divergence  and  vorticity 
fields  since  in  tha  velocity  fields  there  exists  the  sum 
of  sine  and  cosine  terms  while  the  divergence  and  vorticity 
fields  are  represented  by  sine  and  cosines  respectively. 

(In  Table  I,  found  in  the  Appendix,  the  vorticity  and 
divergence  equations  are  listed  term  by  term.) 


As  a first  order  approximation  assume  the  vertical 
velocity,  w,  has  the  same  periodic  distribution  as  the 
negative  divergence.  The  terms  then  may  be  viewed  as 
products  term  by  term  of  a Fourier  series.  If  the  terms 
are  in  phase  then 


/S1 

J o 


in20d0  = 


cos20d0 


and  their  azimuthal  averages  are  non  zero.  If  they  are 
out  of  phase  then 


sin0  cos  9d0  = 0 


and  their  azimuthal  averages  are  zero,  i.e.,  there  exist 
no  positive  and  negative  covariances  from  the  asymmetric 
terms  which  may  contribute  substantially  to  the  hurricane 
budget . 


In  spiral  coordinates  the  integral  over  2tt 
radians  around  the  closed  curve 

of  the  real  or  imaginary  part  of  the  function 


e ( Sg  + i Sr) 


can  be  shown  to  be  zero. 
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may  be  expressed  by 


enS®cos(n  Sr). 
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The  integral  over  the  closed  curve  above  is  given  by 
2tt  2tt 


/ / 


en^®cos(n  S )r  d S r d 

r r 0 


In  polar  coordinates , where  ot  = 0,rdS  rdSA  reverts 

r o 

r d tnr  r d 0,  or  r d r d 9. 

By  substituting 
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r2  = exp  2(Sf  cos  a - SQ  sin  a) 


into  (27)  we  obtain 
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(27) 


to 


(28a) 


(28b) 
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Integrating  SQ  between  0 and  2n  we  obtain 


27r(n-2sin  a)  (n-2sin  a)  f 2S  cosa 
e ^e  e r cos(n  Sr)dSr  (29) 


n-2sina 


Letting  p = n and  a = 2 cosa  and  x = Sr  and  remembering 
that 


ax  , ax,  , 

e cospxdx  - e ( a cos  px  + p sin  px , 

a2  + p2 


we  obtain 


2-r;(n-2  sin  a)  (n-2  sin  a)  2S  cos  a 
e -e  e x 

n-2sin  a 


(2cos  a cos(nS  ) + n sin(n  S ))' 
r r 

4 cos2  a + n2 


The  same  procedure  obtains  for  the  imaginary  part  of  the 
function.  The  importance  of  this  fact  is  that  azimuthally 
integrated  terms  of  a simple  Fourier  series  make  no  con- 
tribution to  the  overall  storm  budget.  Unless  the  term 


t 


is  the  product  of  a Fourier  series  with  the  terms  in  phase 
as  either  the  products  of  sines  or  cosines  as  indicated 
above  it  makes  no  contribution. 

Hence,  if  we  look  at  the  product  of  two  terms, 
such  as  occurs  in  - £6  in  the  vorticity  equation  or  C2 
in  the  divergence  equation  then  we  must  consider  the 
products  of  the  whole  term,  mean  and  Fourier  components  as 
in  (21a)  and  (21b).  The  investigation  of  the  products  of 
terms  A and  B in  the  relevant  equations  gives  us 

(A  + A’)  (B  + B')  = A B + A B'  + A ' B + A ’ B ' 

where  A is  the  azimuthal  average  of  A and  A'  is  the 
Fourier  component.  Both  AB'  and  A'B  vary  in  a simple 
sinusoidal  fashion.  A B does  not  vary  at  all,  and 
A'B'  varies  as  the  product  of  sinusoidal  functions. 

Thus  A B'  and  A'B  contribute  nothing  to  the  azimuthally 
integrated  storm.  5 B,  the  product  of  the  means,  contri- 
butes to  the  symmetric  part  and  A'B'  the  product  of  the 
variances  contributes  to  the  asymmetric  part. 

This  would  at  first  appear  contrary  to  customary 
usage  of  perturbation  theory.  The  underlying  assumptions 
here,  of  course,  are  that  the  asymmetries  are  not 


negligibly  small  in  comparison  to  the  symmetric  part  of 
the  storm. 

The  asymmetric  forms  of  the  terms  for  the 
divergence  and  vorticity  equations  given  in  Table  A1  in 
the  appendix  are  indicated  in  Tables  I and  II  below, 
assuming  that  the  specific  volume  varies  as  the  divergence 
and  the  pressure  varies  as  the  vorticity. 

It  appears  that  as  in  the  symmetric  case, 
advective  and  tilting  terms  may  balance  one  another, 
while  the  convective  and  enstrophy  terms  sum  to  be  equal 
and  opposite  to  the  energy  flux  divergence.  As  opposed 
to  the  symmetric  case,  where  the  brake  on  storm  develop- 
ment appears  to  be  the  frictional  terms,  in  the 
asymmetric  budget,  the  brake  on  the  concentration  of 
kinetic  energy,  the  energy  flux  convergence  arises  out  of 
the  creation  of  divergence  by  enstrophy  and  the  export 
of  convergence  by  the  convective  terms.  The  advective 
and  tilting  terms  appear  uncoupled  from  the  other  terms. 
Possible  subsets  of  equations  in  the  asymmetric  case  of 
the  divergence  equation  are  listed  below. 
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Only  the  advective,  twisting  and  solenoidal  terms  contribute. 
They  must  sum  to  zero  and  be  in  phase  for  a steady  state 
storm. 


The  advective  term  creates  convergence  in  most 

of  the  storm.  The  wind,  U„  .must  be  sufficiently  high  in 

b9 

order  for  the  tilting  term  to  balance  it.  In  the  tropical 

storm  where  Ug  is  trivial,  even  near  the  core,  and 
r 

enstrophy  and  tilting  are  still  at  a low  value  compared 
to  the  mature  storm,  the  convective  term  may  provide  the 
braking  influence  on  the  advective  term.  Since  the  con- 
vective term  is  related  to  the  magnitude  of  the  low  level 
convergence,  a measure  of  storm  intensity,  the  asymmetries 
would  need  to  be  well  developed  for  this  balance  to  occur. 
Otherwise,  the  advective  term  would  lead  the  convective  in 
storm  intensification. 


r 


1 


Again,  let  it  be  stressed  that  this  is  a highly 
hypothetical,  partial  model  of  the  two  dimensional  diver- 
gence and  vorticity  fields  in  the  tropical  cyclone.  Physi- 
cal scenarios  are  suggested  by  a combination  of  phenomenology 
and  the  mathematical  possibility  that  both  the  divergence 
and  vorticity  fields  satisfy  Laplace's  equation  in  both 
symmetric  and  asymmetric  cases.  Thus  far  it  has  only  been 
shown  that  the  vorticity  field  satisfies  the  symmetric 
version  in  a large  majority  of  storms  studied.  Supporting 
evidence  for  asymmetric  satisfaction  by  divergence  and 
vorticity  of  Laplace's  equation  and  symmetric  satisfaction 
by  divergence,  while  indicative,  is  not  yet  firm  enough  to 
warrant  definite  conclusions.  Thus  far,  the  model  is 
phenomenologically  realistic  and  seems  to  be  internally 
consistent  despite  its  elaborate  ramifications. 


The  model  may  be  tested  numerically,  using  the 
model  as  an  analytic  parameterization  of  a fully  developed 
numerical  model  including  energy  and  thermodynamic  consid- 
erations. The  partial  analytical  model  here  deals  only 
with  the  divergence  and  vorticity  equations,  forms  of  the 
equations  of  motion,  and  a simple  form  of  the  continuity 
equation.  The  model  does  not  address  the  first  law  of 
thermodynamics,  the  equation  of  state  or  the  conservation 
) equation  for  the  water  substance.  Further  corroborating 
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evidence  may  be  obtained  by  a Fourier  analysis  of  available 
spiral  data  of  hurricanes  from  satellites  such  as  the 
microwave  detection  of  water  droplet  concentration  in  spiral 
rainbands  studied  by  Allison,  et  al . , (1974).  A study  of 
phenomena,  such  as  water  droplet  content  which  may  be 
correlated  to  the  divergence  field  may  prove  to  be  indicative 
of  the  relevance  of  the  mathematical  solution  to  the  actual 
processes  of  the  storm. 

6.  Conclusion 

X The  generation  of  divergence  and  to  a lesser 
extent  vorticity  are  highly  dependent  upon  the  develop- 
ment of  asymmetries  in  this  model.  Because  of  the  labo- 
rious matnematical  analysis,  a numerical  model  should  be 
parameterized  using  the  asymmetries  of  this  model  to 
assess  their  contributions  to  the  heat,  momentum  and  water 
vapor  budgets  of  the  hurricane. 
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Table  A.l  Comparison  of  Vorticity  and 

Divergence  Equations  by  Terms 


Vorticity  = 

n-  D 6 

Divergence  g—  : 

1) 

Advective  -U*V£ 

Advective  -U*V£  xk 

2) 

Convective  -w 

o Z 

3 6 

Convective  -w  -r— 

3 z 

3) 

aH 

Twisting  -Vw*  — xk 

3U 

Tilting  -Vw*  -r— 

dZ 

9z 

4) 

Divergence  -6c 

Enstrophy  + c2 

5) 

Solenoidal  -VaxVp*k 

Toroidal  -Vct*Vp 
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APPENDIX  A 

A.l  The  Role  of  Terms  in  the  Axially  Symmetric  Vortex 

A. 1.1  Vorticity  and  Divergence  Equations 

In  the  axially  symmetric  vortex,  the  interactive 
terms  of  the  vorticity  and  divergence  equations  come  into 
play  to  reinforce  the  regime  structure.  These  terms  are 
listed  and  named  in  Table  A.l  where  both  equations  are 
presented  vertically  (read  down)  each  term  being  named 
(V  is  the  two  dimensional  gradient  operator). 

The  first  and  third  terms,  which  are  difficult 
to  visualize  because  they  involve  the  cross  product  with 
the  unit  vertical  vector,  k,  are  illustrated  in  Figure  A.l 
for  observed  conditions  under  which  they  are  positive  and 
negative . 
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A. 1.2  The  Characteristic  Equations  for  Given  Regimes 

In  the  vortex  model  ordinarily  there  are  more 
unknowns  than  there  are  equations.  This  entails  guessing 
at  the  values  of  some  of  the  unknowns  in  order  to  obtain  a 
consistent  and  coherent  model.  In  this  section  we  shall 
see  how  we  may  increase  the  number  of  equations  through 
grouping  analytically  similar  terms  which  should  therefore 
sum  up  in  a self  consistent  way  to  satisfy  the  model. 

Since  the  vortex  is  modeled  in  terms  of  vorticity 
and  divergence  distribution  the  model  gives  an  indication 
of  how  the  individual  terms  in  the  vorticity  and  divergence 
equations  behave.  For  the  steady  state  vortex  the  terms 
in  scalar  form  appear  to  fall  into  discrete  analytical 
categories.  The  categorization  is  based  on  the  terms  having 
similar  analytical  functions  of  radius  and  are  of  three 
types.  These  are  given  in  Table  A. 2.  The  categorization 
is  a function  of  regime  and  is  given  in  Table  A. 3 for  the 
three  regimes. 

In  Table  A. 3 the  terms  are  listed  by  class  for 
the  vorticity  and  divergence  equations.  The  sign  of  the 
term  is  given  in  parenthesis.  In  the  core  regime,  the 
Class  II  terms  in  the  vorticity  equation  change  sign  from 
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the  "eye"  portion  of  the  core  to  the  "eye  wall"  and  hence 
have  a sign  for  each. 

Assuming  the  coefficients  of  Class  I sum  to  Af(r), 
Class  II  sum  to  Bg(r)  and  Class  III  sum  to  Ch(r)  then  the 
steady  state  equation  may  be  written: 

Af(r)  + Bg(r)  + Ch(r)  = 0 (A.l) 

In  the  core  C is  zero  since  V2s  = -VaxVp  = 0 so  that 

2 

Af ( r ) + Bg( r ) = 0 (A. 2) 

If  A and  B are  non  zero  then  Af(r)  = Bg(r).  But  this 
means  that  f(r)  - -CQg(r)  where  CQ  is  some  other  constant. 
This  may  be  true  at  one  point  but  is  not  true  elsewhere. 

Thus,  if  A and  B are  zero  f(r)  f CQg(r)  since  CQ  = -B/A 
is  undefined. 

Assuming  the  parallel  terms  in  the  divergence 
equations  behave  in  the  same  fashion,  then  the  terms 
-Va*Vp  - V2K  must  be  of  the  form  Ch(r)  = 


0. 
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Therefore  there  exists  a whole  series  of  term 
couplings  which  shed  light  upon  the  vertical  distribution 
of  the  terms.  Thus,  there  exist  the  equations: 


(-) 

-U*7? 


Vw-  xk  + -gi-  v = 0 

9z  9z  z 3z 


(A. 3) 


-w  It  - - 0 


(A. 4) 


Since  we  may  solve  (A. 4)  to  obtain  = c6/w 

oZ 


and  we  know  that  UQ  = — ^ f £ dr  so  that 

tt  «s7T  J 


3U0  1 r r 9t 

■jzj-  = 2^  / 9z  then  we  obtain  by  substituting  (A. 4) 


into  (A. 3) 
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(A. 5) 


(A. 6) 


We  therefore  may  obtain  the  value  of  the  vertical  coefficient 
of  eddy  viscosity  in  the  core. 


By  specifying  the  analytical  form  of  the  scalar 
terms  for  a given  regime  we  may  solve  for  variables 
otherwise  unattainable  such  as  the  eddy  viscosity,  v , as 
a function  of  height  and  the  specific  volume,  a,  as  a 
function  of  radius.  In  this  way  we  are  able  to  diminish 
the  number  of  unknowns  in  the  model  because  we  have  increased 
the  number  of  equations. 
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Appendix  B 

Creation  of  Regimes  in  the  Vortex 

B.l  Physical  Parameters  Governing  the  Vortex  Structure 
Perhaps  the  easiest  way  to  understand  the 
rationale  for  the  lateral  distribution  of  the  regimes  of 
the  vortex  model  is  to  follow  the  suggested  sequence  of 
events  which  lead  to  the  formation  of  the  vortex. 

In  order  to  do  this  it  is  necessary  to  have  an 
understanding  of  the  role  of  the  vorticity  and  divergence 
equations  in  the  makeup  of  the  model.  Moreover,  the  role 
of  the  Richardson  number  is  essential  to  understanding 
how  order  may  be  imparted  to  an  otherwise  amorphous 
structure . 

B.2  Richardson  Number 

The  Richardson  number  expresses  the  characteristic 
ratio  of  the  energy  of  a work  done  against  gravitational 
stability  to  energy  transferred  from  mean  to  turbulent 
motion.  Theoretical  studies  have  placed  the  critical 
Richardson  number  from  .25  to  2 with  stability  for  higher 
values  and  instability  for  lower.  The  Richardson  number 
may  be  expressed  by 
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atmospheric  oceanic 

Since  air  spiralling  into  the  hurricane  must 
cross  isobars  to  lower  pressure  the  air  experiences 
cooling.  Air  within  the  cloud  layer,  however,  cools 
moist  adiabatically , i.e.,  more  slowly  due  to  the  release 
of  the  latent  heat  of  condensation  of  water  vapor  while 
air  below  the  cloud  cools  dry  adiabatically.  Since  water 
vapor  liberates  597.3  calories  per  gram  upon  condensation, 
the  lower  air  cools  much  more  rapidly  than  the  upper.  The 
brake  in  the  process,  if  the  vortex  is  a hurricane,  occurs 
due  to  transfer  of  sensible  heat  from  the  sea  surface 
into  the  vortex  boundary  layer.  The  resulting  stability 
is  only  counteracted  by  increased  turbulence  due  to  a 
higher  mean  velocity  as  one  approaches  the  center  of  the 
vortex.  The  more  vigorous  mixing  tends  to  destroy  the 
inversion  caused  by  the  differential  cooling.  It  is  thus, 
however,  that  the  vertical  velocity  field  tends  to  be 
concentrated  into  a narrow  throat. 
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The  moist  expansion  in  the  cloud  layer  leads 
to  a higher  potential  temperature,  0 of  the  air  than  its 
lower  neighbor.  This  is  expressed  by 

0g  = 0 exp 

where  L is  the  latent  heat  of  evaporation,  597.3  cals/gm, 

w the  saturation  mixing  ratio  of  the  water  vapor,  C 
s p 

the  specific  heat  of  dry  air,  and  T the  temperature  of  dry 
air . 

B.3  Sequence  of  Events  in  the  Formation  of  a Severe 
Atmospheric  Vortex 

Given  an  arbitrary  distribution  of  vorticity, 
divergence  and  pressure,  the  following  sequence  of 
events  is  hypothesized  as  leading  to  the  formation  of  a 
severe  vortex. 

Assume  an  initial  parabolic  distribution  of  pressure, 
divergence  and  vorticity,  with  a subsequent  radial  and 
tangential  velocity  distribution  as  shown  in  Figures 
Bla  and  Bib.  Assume  also  that  the  total  amount  of  vor- 
ticity and  divergence  are  constant  but  may  be  and 
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indeed  are  rearranged  in  response  to  the  influence  of  the 
various  terms  of  the  relevant  equations. 


V 
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The  first  step  in  this  process  is  the  rearrange- 
ment of  the  distribution  by  the  two  dimensional  Laplacians 
in  the  vorticity,  divergence  and  pressure  tendency  equations. 
The  distributions  then  become  as  illustrated  in  Figures 
B2a  and  B2b . 
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The  second  modification  occurs  when,  in  response 
to  the  concentration  of  the  wind  peaks  inward,  the 
inspiralling  air  experiences  differential  cooling.  The 
extent  of  the  cooling  depends  significantly  upon  whether 
the  trajectory  of  the  inspiralling  air  is  within  a cloud 
or  below  it.  This  differential  cooling  causes  an  inversion 
to  form  faster  than  turbulence  can  break  it  up  in  the  outer 
regions  of  the  vortex.  This  in  turn  suppresses  vertical 
velocity  making  the  formatiom  of  the  inversion  even  more 
pronounced.  Thus  the  critical  Richardson  number  Ri 
defines  the  boundary  between  the  inner  and  outer  regimes 
(Figures  B3a,  B3b). 

The  final  stage  occurs  when  the  enstrophy  term, 

C2,  becomes  so  large  at  small  values  of  radius  (due  to 
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crowding  of  the  vorticity  and  divergence  into  a smaller 
and  smaller  area)  that  divergence  is  created  so  fast  by 
this  term  that  the  vorticity  peak  collapses  like  a 
volcanic  caldera  bringing  about  the  final  stage  in  the 
cycle  illustrated  in  Figures  B4a  and  B4b  with  the  creation 
of  a core  interior  to  the  inner  regime. 

The  fit  of  the  model  to  the  hurricane's  pressure 
profile  is  also  very  good  and  substantiates  the  presence 
of  the  three  regimes  mentioned  above.  An  extension  of 
the  concept  which  accounts  for  the  pressure  and 
tangential  velocity  fields  presents  a logical  explanation 
for  the  radial  and  vertical  velocity  fields. 


» 


The  organization  of  the  storm  from  a preexistent 
but  amorphous  field  of  vorticity  and  convergence  into  the 
highly  structured  vortex  is  seen  as  the  effect  of  turbulent 
lateral  diffusion  constrained  by  critical  values  of  the 
Richardson  number  and  of  enstrophy  (the  square  of  the 
vorticity).  The  interaction  of  these  mechanisms  provides 
a plausible  scenario  for  the  structuring  of  the  mature 
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storm. 

changing 

insights 


The  application  of  these  concepts  to  the  ocean 
what  must  be  changed will  provide  valuable 
into  the  behavior  of  the  upper  ocean. 
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Figure  la  and  lb.  Vorticity,  divergence 
and  velocity  regimes  in  the  zero  Laplacian 

vortex . 
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Hurricane  Cleo 


Figure  2a.  Piecewise  continuous  fit 
of  Zero  Laplacian  Vortex  to  Riehl’s 
(1963)  hurricane  Cleo  data.  Rankine's 
combined  vortex  is  included  for  comparison 
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Figure  2b.  Piecewise  continuous  fit  of 
Zero  Laplacian  Vortex  to  Riehl's  (1963) 
hurricane  Hannah  data.  Rankine ' s combined 
vortex  is  included  for  comparison. 


Figure  3.  Rotation  of  Cartesian  Coordinate 
System  through  angle  a. 


Eye  Wall  Hurricane  bar 

Figure  4.  Rotation  of  semilog  coordinate  system 
Inr,  9 to  produce  log  spiral  system 


Spiral  Space  with  angle  of  rotation  20°. 
The  spirals  occur  at  every  - y . The 
inner  values  of  S„  would  be  completely 
eclipsea  by  a strong  axially  symmetric 
phenomenon  at  the  storm  center. 


